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I. INTRODUCTION 

In practice, the control system designer rarely has the 
freedom to feed back the entire state of the system to 
achieve the desired system performance. Thus, an important 
problem from a practical standpoint is the determination of 
constant, output feed back gains for the control of systems 
with unmeasured states. 

The pole placement only and pole placement and robust- 
Less design procedure uses placement to establish a designer 
selected performance level and then a minimum singular value 
level to establish robustness. The pole placenent only was 
relatively easy to implement through a numerical optiniza- 
tion routine. By using this numerical procedure it is also 
Simple to incorporate robustness into the procedure along 
with the performance requirement. The technique, which 
utilizes a modern optimization routine, can significantly 
assist the designer in obtaining robustness in the face of 
Cross coupling perturbations. It has been shown that the 
cross coupling problem can be detected by using classical 
open loop Bode diagrams as well as modern control analysis 
based on singular values. As currently employed, the pole 
placement and robustness program is used to obtain pole 
placement and robustness for a given set of starting gains 
and a selected optimization routine from the ADS progran. 

The pole placement and robustness design procedure is a 
straight forward numerical optimization procedure for the 
practical application modern MIMO system analysis. The new 
aspects of the procedure are the implementation of both pole 
placement and robustness criteria within the same design 


program. 
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The pole placement and robustness design routine devel- 
oped for this thesis has been used on X-22A A/C problem. In 
these studies the pole placement and robustness design code 
125 proven capable of meeting the desired goals of pole 
placement and robustness and also brought to light some 
interesting aspects of the cross-coupling perturbation 
problen. 

The remainder of the thesis will present background 
material on Robustness (SISO system) in Chapter Two and 
Robustness of MIMO systems in chapter Three. Pole piacement 
only design procedure in chapter Four, Pole placement and 
Robustness design in chapter Five, along with a discussion 
of Applications for X-22A a/c problem by input transfer 
function F*G analysis and output transfer function G*F anal- 
ysis. Conclusions will be presented in the final chapter. 
The computer programs used in the present analysis vere 
developed in (Ref. 9). 
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II. ROBUSTNESS: SINGLE-INPUT SINGLE-OUTPUT SYSTEMS (SISO) 

A review of the concept of robustness and stability in 
the framework of a conventional SISO system will be done 
before pursuing the concepts ina more complicated fashion 
in the following chapters. A simple interpretation of 
robustness is the ability of the system to tolerate design 
perturbations. These perturbations could be in the form of 
actuator failures, plant parameter uncertainty, unmodeled 
dynamics or nonlinear terns, or any one of many other 
perturbations to the nominal design of the system. 

The primary reason for feedback systems is the control 
of uncertainty within the system. By appropriate use of 
feedback, properties that would lead to an unstable system 
may be controlled. When stability and robustness aspects 
are considered for a SISO system, frequency iomain design 
concepts, using either Nyquist or Bode plots, are normally 
used. Robustness in SISO systems is formulated naturally by 
tLe concept of gain and phase margins, both of which are 


readily available on the Nyquist or Bode diagram. 
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Figure 2.1. Classical Bode Plot. 


can be defined in terms of the 


in either the Bode or 


Gain and phase margin 


open-loop frequency domain 


Nyquist format. Figure 2.1 depicts a classical Bode plot 


showing gain and phase margin determination from the plot. 


used to obtain this informa- 


that if the open-loop 
poles in the right halt 


The Nyquist plot may also be 


+ Nyquist criterion 


transfer function  G(s)H(s) 


plane and the limit of G(s)H(s)-constant as s—oothen for a 
stable system the locus of G(s) H(s) Will encircle the -1+j0 
point n times in the counterclockwise direction aS s varies 
along the Nyquist contour. If there are no poles in the 
right half s plane then the locus will not encircle the -1 + 
jO point. The diagram in Figure 2.2 illustrates a nominally 
stable system. The gain and phase margin may be determined 
directly from the diagram. 

Any change in the loop transfer function, provided the 
order of G(s)H(s) does not change, that changes the number 


of times the locus of G(s)H(s) encircles the (-1,0) point in 
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Figure 2.2 Nyguist Plot of Stable Systen. 


the Nyquist plot causes the system to become unstable. This 
leads to the conclusion that the minimum distance of the 
locus of Gí(syH(s) to the (-1,0) point is a measure Of the 
system stability. This distance concept carries over 
directly to the Multi-Input Muiti-Output (MIM3) system as 
will be shown in the next chapter. Examples of a multipli- 
cative perturbation and an additive perturbation illustrate 
this idea. Pigure 2.3 is an additively perturbed systen. 
Figure 2.4 shows the Nyquist plot for this system. Assuming 
that the plant is itself stable and the perturbations are 


also stable the diagram may then be used to determine how 
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Figure 2.3 Additively Perturbed Systen. 


near the system ís to instability for the given perturla- 
tron. 

Since the system is stable the (-1,0) point is encircled 
the correct number of times by the nominal plant. If the 
locus of g(jW) in the diagram is warped until it passes 
beyond the (-1,0) point then clearly the number of encircle- 
ments of this point will change and the system will become 
unstable, assuming the order of the plant is not changed by 
the perturbation. To keep the locus of points from moving 
beyond the (-1,0) point equation 2.1 must hold. 
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Figure 2.4 Additive Nyquist Plot. 


LAgG(Jw) I<! 1+ g (Jw) I] (2.1) 


This condition is illustrated in Figure 2.5. The right-hand 
Side of equation 2.1 is just the magnitude of the return 
difference transfer function of the nominal systen. The 
multiplicative case is depicted in Figure 2.6 with its asso- 
ciated Nyquist plot in Figure 2.7. The recuirement for 
Stability is similar to the additive case and may be stated 


in equation 2.2 


ıAg(jw)I <| 1 + (g(jw) 71) ] (2.2) 
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Figure 2.5 Nyquist for Inequality Additive Condition. 


The above arguments will be applied again in Chapter 3 to 
develop multivariable stability and robustness properties. 
With this basic review of the concepts of stability and 
robustness in the classical SISO system complete, the next 
Chapter will extend some of these basic concepts to the AIMO 


System. 
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Multiplicative Systen. 





Figure 2.6 


—————————É— ا‎ GUÐ Sl SNIÐ A n" — سس‎ A — 


-1 pe 6 
^ 

S 
کس سے کا‎ 


0 ۲ 


E 


Figure 2.7 Nyquist Plot for Multiplicative System. 


III. MULTIVABIABLE SYSTEMS ROBUSTNESS 
A generalization of the SISO Nyquist theory discussed in 
the previous chapter has been made for the MIMO proLblem. 
This generalization leads directly to the apjlicatior of 
Singular value concept. The generalization is expressed in 
the form of the multivariable Nyquist theorem which requires 
that a closed loop stable system have the same number of 
counterclockwise encirclements of the origin by the locus of 
the det (I+G(jw)) as the number of open loop poles that are 
in the right half plane. This theorem is formally stated as; 
Let N[f(s)] denote the number of clockwise encirclements of 
(-1,0) by the locus of f(s) as s traverses the contour D 
of Figure 3.1 in a clockwise sense. The closed-loop system 
will Le stable if and only if for all R sufficiently large 
N[f(s)] » -P 


where N = number of encirclements 
f(s) - -1 * det( I*G (s) ] = 1, 이. 1 Tae 
P = the number of closed right-half plane zeros of 


7 7 


The application of the Nyquist theorem comes through the 
fact that a multivariable system will not be robust to 
modelling errors if the return difference matrix, I + G, is 
nearly singular for some frequency. If I + G is nearly 
singular a small change in G may make I + G exactly 
Singular. This causes the det(I + G) to become zero and the 
Nyquist encirclement count to change indicating an unstable 
system. It is possible for very small changes in ItG to 
produce large changes in the determinantof I+G. The matrix 
1+0 
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Figure 3.1 Nyguist D Contour. 


10/s*a 9+95 
10/sta 10/sta 
has determinant 0.1/(sta)?. If the element p 1s changed 


by only one percent to 9.9/sta the determinant becomes 
1.1/ (sta)? which is a significant change in the determinant- 
value. Therefore, it is evident that det (Il + G) is not an 
accurate measure of how near the return difference is to 
Singularity. Researchers, in the field of controls 
[Re£. 1], [BRef. 2], (Ref. 3], (Ref. 4] have used singular 
value analysis to determine how near the return difference 
matrix is to singularity. 

Since the number of encirclements of the Nyguist diagran 
changes as f(s) passes through the -1 point or when det (I+G) 
is zero it is important to find how near the return differ- 
ence matrix I+G is to being singular. This nearness to 


Singularity can be interpreted as the distance of matrix I*G 
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to the critical point, -1. A Quantity woteh can be used 
express the nearness to singularity of the matrix is the 
minimum matrix singular value denoted by Q^ . Given a matrix 
A the singular value may be expressed by equation 3.1 


en ë 


O\ (A) = min ( JN, (A A)) (3.1) 
4 


where A: (dE) is the eigenvalue of the complex conjugate 


transpose of A times A. A kasic MIMO linear system is 


ke eee 


nF EE a Ere oo eS eee eee 


Figure 3.2 Basic Multi-input Multi-output System. 


depicted in Figure 3.2. An additive perturbation to the 
plant is shown in Figure 3.3. If the plant is stable before 
the perturbation is added to the system the Nyquist theorem 
will be satisfied and the locus of GH will not encircle the 
-T70 erreal pont - When the perturbation is added to the 


system as long as the Nyquist locus is not forced to 
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Figure 3.3 Additive Perturbation. 


encircle the -1 point the system, will remain stable A 
sufficient condition, recalling the SISO discussion in 
chapter 2, for the perturbed Nyquist plot not to change 
encirclements is that the norm of the perturbation AG 
remain less than the norm of the return difference matrix as 
expressed in equation 3.2. 


11 AS(w)II« 1711 (1*8) =) || (3-2) 
The condition W 20 will that the locus of the det (I*G) 
does not pass through the -1 point. If the 1 or Euclidean 
norm is assumed for this condition the equation 3.2 may be 


expressed in terns of singular values as equation 3.3. 
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— ss 


OW (AGS ® (I + €) (3.3) 


This result states that as long as the maximum singular 
value of the perturbation matrix ÁG is below the minimum 
norn value of the return difference matrix the system will 
remain stable. The problem of guaranteeing robustness 
becomes that of finding the largest norm of the perturbation 
quantity, the largest singular value, for whicn the smallest 
norm or singular value of the return difference matrix will 
remain non-singular. 

The multiplicative form for a system such as Figure 3.4 


gives the similar norm equation in eguation 3.4 


I1 G (ju) IIX1/101 (I + (S) 7) 11 (3.4) 


0) 20. which may be expressed as 
QUAS $ QUIS G7) (319) 
Singular value decomposition software 1s readily available 


to determine how near the matrix ItG or I+ (6) is to 


Singularity. 
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Figure 3.4 Multiplicative Perturbation. 


IV. POLE PLACEMENT ONLY DESIGN PROCEDURE 

The designer may uSe either objective, constraint or a 
combination of functions to secure the desired pole loca- 
tions. As currently implemented in the program the cost or 
objective portion of the pole placement procedure is 


constructed as equation 4.1 


£ 
OBJ =2 (68 Ar)? + ۱ ٣ (4.1) 
i71 ۱ | 
where ^, = real eigenvalue 
Ke imaginary eigenvalue 
Ng = desired eigenvalue location 
Ni = desired eigenvalue location 


The constraint fcrmulation is a function that must be 
kept negative or the constraint is violated. It is written 


aS equation 4.2 


O(N = (Nr ۸" ) 4 ٣ (4.2) 


where r is a tolerance circle established as a function of 
pole placement position. Since the aim of the optimizer is 
to keep g negative any time the N function of the 
constraint is greater than r the constraint will become 
active, i.e. violated. The optinizer will then attempt to 
nove the constraint to the inactive status by adjusting the 
design parameters of the system. For obtaining the feedback 
gains and desired poles use the 'CONXSV' program (Ref. 9], 
weight function 1 equal 1 the other weight functions are set 
to 0. After obtaining the feedback gains the NPS 'OPTSYIS' 
program was used to plot the Bode, Nyquist diajram and pole 


zeros plot. 
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Ve POLE PLACEMENT AHD ROBUSTNESS DESIGN PROCEDURE 

Consideration of implementation of the freguency domain 
or robustness portion of the design procedure begins with 
the concept of MIMO phase and gain margins. Several useful 
theorems on Singular value analysis of multiloop systems are 
presented in (Ref. 3]. One of these theorems relates the 
matrix Singular value of the return difference function toa 
parameter,x , and further shows that as long as the maximum 
singular value of the perturbation function (L-1 - I) 
remains less than this, the system renains stable. The 
value oi d is then related to gain and phase margins of the 
MIMO system. The relationship developed is given in ejua- 
Mons 5.1 and 5.2: 


gain margin = GM =1/(1+0,) )5.1( 
phase margin = PM = *cos-!(1- X2/2) (5:2) 
provided that equation 5.3 holds. 


Q(I*G)2 de [5* 3) 
for some do <1 
These phase and gain margins are guaranteed in every loop 
simultaneously. 
Universal gain and phase margin curves, (Ref. 5], based 
on the minimum singular values of the return difference 


Matrix are developed from equation 5.4. 


(L-1-I)-max d 뼈 - 121, ) =, (1-205 0) (5.4) 
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for all n with 0. These curves shown in Figure 5.1 
allow the designer to pick a singular value that corresponds 


to a specific gain and phase margin for a given system. 
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Figure 5.1 Universal Gain and Phase Singular Value Plot. 


Since the universal curve in Figure 5.1 provides a 
convenient method of specifying gain and phase margins in 
terms of singular values the robustness portion of the pole 
placemert and robustness design procedure uses the minigun 
Singular value level of the return difference matrix to 
determine the robustness. The minimum Singular value level 
is formulated as an objective or constraint function in 


equation 5.5 


J= ) (max (0, tÜ5 - c (ju, p)))? (5885) 
A 
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The optimization procedure may be used to change feedback 
gains until the minimum singular value ís raised above this 
desired design level. Although the same formulation can be 
used as a negative constraint function it has not been 
implemented as such within this program. There are numerous 
ways the singular value formulation could be implemented 
within the program by changes of the code if design require- 
ments forced such changes. 

The pole placement and robustness design program is 
hased onthe ADS code to implement the design variable 
selection procedures. The pole placement and robustness 
program is used to provide designs for state or output feed- 
back problem. 

For the state or output feedback design program the user 
nust input the plant matrices A,B,C and initial starting 
Mates tor the feedback matrix F. The matrices correspond to 


the following linear differential system: 


x = Ax + Bu (5.6) 
y = Cx (5.7) 
u = -FX (5.8) 


As the design program is currently coded the user may 
run output feedback or state feedback design by specifying 
the C matrix as the diagonal(I) matrix for state feedback. 
The program relies on initial starting values of the feed- 
Pack gains, E. 

The ability to select acceptable starting values for the 
feedback gains will make the procedure more efficient in 
operation. As currently employed, the program is used to 
obtain pole placement and robustness for a given set of 


Starting gains and a selected optimization routine from the 
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ADS program. If the optimizer is not able to meet the 
desired design goals on this program run two options are 
available. First, change to a different optimization routine 
from the list of available ADS routines and rerun the 
problen. This was usually successful in improving the 
design. Second, the designer useS a new set of starting 
values for the feedback gains and repeats the design proce- 
dure. Both options might be used on particularly difficult 
Cases. 

The pole placement and robustness design algorithm 
computes input additive, output additive, input multiplica- 
tive, and output multiplicative singular values. The versa- 
tility of the pole placemeit and robustness design is 
obtained by incorporating a state of the art optimizer 
routine ADS, with currently available singular value compu- 


tation routines. 
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VI. APPLICATIONS(X22A LONGITUDINAL PROBLEM) 
This chapter will deal with a more practical application 
of tne numerical optimization program. In this problem the 
combined pole placement only, robustness design procedure 
wiil be applied to the linear longitudinal dynamic channels 
of an X22A V/STOL a/c [Ref. 8]. The model is that of the 
longitudinal dynamics of an X22A V/STOL a/c at low altitude 
and airspeed = 65 knots. The dyrnauic model of tne svstem is 


E-CAx*5 (6.1) 


ic 


X s(u,u,q,8) | (6.2) 
u = (rd) (6.3) 
where 
#0. 15 003 ا‎ 7 
A = EZ =02 59 109. 43 2298 
-0.01 - 2 = Dos, 0 
E. 0 0 le 0 
and where 
209305 0.52 
B = 0% =], 
035 0.02 
0. 0. 


with full state available for feedback. Table 1 is a summary 
of parameters. The system is not open-loop stable. 

The control law is formulated to  satisry the desired 
performance specification. Equation 6.4 is the basic control 


law. 
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TABLE 1 
X22A V/STOL a/c Parameter Definitions 


Body axis forward velocity perturbation 













ft/sec 
ft/sec 


rad/sec 





Vertical velocity perturbation 









Pitch rate 
Pitch attitude 


Elevator position 






rad 






inch 






inch Throttle position 


(6.4) 


~ 


——————-— m مہ‎ 


Variable 


design studies presented up to this point have been 
based on breaking the system locp at the input as shown in 


2 


Figure 6.1 System Block Diagran. 
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Tae 





2200267762 1 In multivariable tueory the location of the 
break in the loop changes the return difference OL tne 
systen and the transfer function formulation. Ir Figure 6.1, 
number 1 depicts a system with an input loop break point 
while number 2 depicts an output loop break point for output 
return difference determination. The return difference func- 
tion for point 1 is written as 1+5 while the return differ- 
ence for point 2 is I+tGF. In figure number 1 the transfer 
function FG is 
FG = F'* (SI - A) 1*B 


FY = PRC 
and GF is GF - C*(SI - A) !*B' 
B' - BX*F 


A.  INPUT TRANSFEB FUNCTION F*G ANALYSIS 


Using the NPGS 'CONXSV' progran [Ref. 9], to get the 
desired feedback constants, the pole placement only progran 


was first run to get the desired pole and tne computed 


TABLE 2 
Pole Only Desired Pole and Computed Pole 


Desired pole Computed pole 
2 





pole(see Table 2). The feedback matrix F is 


B - ZU ا ا‎ 0 1 13. 97235 7.73604 
UR U0. 10997 =-90.41902 12742874 
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Figure 6.2 Pole Placement Only Singular Value Plots. 


the objective function-is 0.0002 and the adiitive minimum 
singular value vs frequency is seen in Figure 6.2. This 
plot shows that pole placement only design has very low 
minimum singular values for the return difference matrix. 
Pole placement only goes as low as -22 db or 0.0827 rad rear 
2 rad/sec in frequency. Using the universal gain and phase 
diagram as discussed in Chapter 5 this equates to a gain 
margin of about 0.95 db to 1.15 db and a phase margin of 4.7 
degrees. ^ These phase and gain margins are quite small, 
Showing the need to run a pole placement and robustness 
design program. It is assumed that the eigenvalues (pole 


locations) of the system, as developed in (Ref. 7], are the 
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required poles for the performance criteria. Once the pole 
locations are set, robustness criteria must be selected. 
From the universal gain and phase Margin curve discussed 
earlier, Singular value levels were selected for this 
problem. The Singular value level chosen was 0.6 rad. This 
corresponds to a gain maryin of -4.0 db to 8 db and a phese 
Margin of about 35 degrees. 

For a singular value level cf 0.6 rad the pole placement 
and robustness design routine places the poles as shown in 


Table 3. The slight differences in these pole locations 


TABLE 3 
Pole Placement and Robustness Design Pole Locations 





have an insignificant effect on the performance. 


The feedback constant F is 


F = -0.07493 -0.08746 8.0529 16222076 3 
OSS 2. 49 3 7.6 4 


The ' objective function is 0.0000 and the additive 
Minimum singular value vs frequency is seen in Figure 0.3. 
This piot shows that in the robustness design the gain 
adjustment moves the minimum Singular value from about 0.08 
rad, with very poor phase and gain margins, toa level of 
0.645 rad, above the desired values of gain and phase. 
After obtaining the feedback gains the NPGS 'OPTSYS' program 
was used to obtain the necessary data, time response plot, 
Bode plot, and pole-zero map. The pole placement only 


design closed loop time response plot is shown in Figure 
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Figure 6.3 Robustness Design Singular Value Plots. 

6.4. It shows good pitch rate and pitch attitude responses 
to de step input, but slow responses to the forward 
velocity perturbation and vertical velocity perturbation to 
d Step input. The input transfer function F*G robustness 
design time response plot is shown in Figure 6.5. This 


indicates a better time response for pitch rate and pitch 
attitude inputs than pole placement only design, also good 
vertical velocity perturbation response to g step input, 
but sluggish response for the forward velocity perturbation 
[Bef. 8]. Comparing the loop Pode plots of the poie place- 
ment and robustness design transfer function, F*G, shows an 


increase in robustness. 
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Pole Placement Only Design Time Response. 
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Figure 6.5 F*G Robustness Design Time Response. 
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Considering a SISC Bode analysis, the channel input 1 
output 1 bode plot are shown in Figure 6.6 and 6.7. 
The bandwidth(BW) shows a small decrease, 3 rad/sec, and an 
increase, in the phase margin(EM), of 22 degrees. he DC 
Sain 02965 fron 34.3 daD to 25.7 db, ^a 9 db decrease and the 
Slope changes very atruptly inthe pole oniy design (see 
Figure 6.4). The robustness design seen in figure 6.5 urops 
the curve down 20 db due to the two zeros at 0.3 rad/sec, 
which move toward the two poles at 1.27 rad/sec, and the 
other zero moves closer to another pole 0.16 rad/sec during 
the optimization. This channel only shows a smoothing of the 
curve. Consider the channel input 1 output 2 bode plots 
Shownir Figures 6.8 and 6.9. This channel shows a large 
shift in the bandwidth (PET) from 30 rad/sec to 2.83 rad/sec. 
The DC gains are similar and increase the gain margin 24 db, 
but decrease phase margin 90 degrees. Figure 6.3 shows that 
the frequency response curve slope changes detrimentally, 
but, the robustness design, as shown in figure 6.9 shows a 
decrease in the slope down 20 db curve by zeros move the one 
pole 0.17 rad/sec for pole-zero cancel andthe other one 
zero moves on the minimum singular value frequency 1 
rad/sec. The reduction of bandwidth and increase in gain 
margin yields an increased tolerance to perturbation. The 
channel input 2 output 1 bode plots are shown in Figures 
6.10 and 6.11. The DC gain decreases 10 db and bandwidth 
shifts left 0.1 rad/sec, increasing the phase margin 43.7 
degrees. The zero shift has the effect of smoothiny the 
frequency response curve in the vicinity of the frequency of 
the minimum Singular value. The channel input 2 output 2 
bode plots are shown in Figures 6.12 and 0.13. These 
figures show the same response as the channel 1;2 but a 
bandwidth shift right of 0.47 rad/sec and the same DC gainas 
for a smooth curve, caused by the zero location shift in the 


sane manner as the other channels(channels 1;1,2;1,2;2). 
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Figure 6.6 Pole Placement Only Design F*G 1:1. 
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Figure 6.7 Robustness Design F*G 1:1. 
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Figure 6.8 Pole Placement Only Design F*G 1:2. 
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Figure 6.9 Robustness Design F*G 1:2. 
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Figure 6.10 Pole Placement Only Design F*G 2:1. 
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Figure 6.11 Robustness Design  F*G 2:1. 
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Pole Placement Only Design F*G 2:2. 
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Figure 6. 12 
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Figure 6.13 Robustness Design  F*G 2:2. 
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The F*G transfer functions bode plot numerical results are 


skown in Table 4. 


TABLE 4 
F*G Transfer Functions Bode Plot Results 


Bandwidth 
rad/sec 
SV 


pole 





The most noticed change in the overall system, however, 
is in the transfer function input 1 output 2. This is the 
change at the minimum singular value position, which greatly 
reduces the gain and bandwidth in this channel through a 
change in feedback gain. The cptimizer routine brings the 
entire system gains to more balanced conditions and recovers 
a highly robust design. 

The yain changes associated with the robustness improve- 
ment cause the zeros of the various closed-loop pole-zero 
diagram of the closed-loop transfer matrices to move. A 
comparison of the eight pole-zero diagrams is shown in 
Figures 6.14 to 6.17. 

The significant feature of these pole-zero diagrams is 
the shift in the zeros of the optimized design in a direc- 
tion that attempts to equalize or balance the frequency 
response for frequencies in the vicinity of the minimum 
Singular values. The pole-zero diagram of channel input 1 
output 2 will be discussed as an example of this effect. In 
figure 6.12 the pole only design zeros are located at about 
-0.31 and -0.78. When the pole placement and robustness 


routine has completed the feedkack gain modification these 
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3 and 2:4. 


Pole-Zero Map for 2 


Figure 6.17 
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zeros have shifted to -0.17 and -0.56. The effect of these 
Zero shifts is tO Combine with the pole locations to 
equalize the frequency response and increase DC gain of 14 
db for the same bandwidth, at the minimum singular value 
frequency 1 rad/sec as depicted in Figures 6.18 and 6.19. 
zero shifts for the remainder of the transfer functions 
provide Similar results in the other channels. By moving 
toward the frequencies associated with the minimum singular 
values the zeros have balanced the overall ‘frequency 
response of the system in eacn channel. While the channel 
gain modification is the primary mechanism for robustness 
recovery, the zero shift associated with the feedback gain 
Cnanges is directly related to the overall frequency 
response of the system. A rokustness design meeting the 
required pole locations and a robustness singular value 
level of 0.6 rad provides adequate gain and phase margin for 


the design. 


B. OUTPUT TRANSFER FUNCTION GFF ANALYSIS 


Using the NPGS 'CONXSV' program [Ref. 9], to get the 
desired output transfer functicn feedback constants, the 
pole placement only program was first run to get the desired 
pole and the computed pole(see Table 2). the feedback 
matrix F is 


F 


053. 0.50675 13.97235 7.73604 
Umu -0.10997 -90.41902 12.42874 


the objective function 15 0. 2 and the additive output 
minimum singular value vs frequency is seen in Figure 6.20. 
This plot shows that pole placement only G*F design has very 


low minimum singular values for the return difference 
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Figure 6.18 Pole Placement Only Design F*G Closedloop 1:2. 
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Figure 6.19 Robustness Design F*G Closedloop 1:2. 
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Figure 6.20 Pole placement only S*F Singular Value Plots. 


matrix. Pole placement only goes as low as -48 db or 0.004 
rad near 1 rad/sec in frequency. Using the universal gain 
and phase diagram as discussed in Chapter 5 this equates to 
again margin of about -0.003 db to 0.003 db and a phase 
margin of 0.5 degrees. These phase and gain margins are 
quite small, showing the need to run a pole placement and 
robustness design program. it is assumed that the eigenva- 
lues (pole locations) of the system, as developed in 
[Ref .2 75% are the required poles for the performance 
criteria. Once the pole locations are set(see Table 5), 
robustness criteria curve discussed earlier, singular value 
levels are made for this problem. The singular value level 
chosen is 0.6 rad. This corresponds to a gain margin of -4.0 


db to 8 db and a phase margin of about 35 degrees. 
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For a singular value level of 0.6 rad the pole placement 
and robustness design routine places the poies as shown in 
Table 5. There isa trade off between performance and 
robustness in many control problems. ۲ (٣٢0.7۴ 
function G*F singular value calculation, the shift in this 
pole can have a Significant effect on the performance, it 
was found necessary to nove the pole at -0.5 to -0.05 in 
order to achieve the desired rokustness. The slight differ- 


ences in these pole locations an insignificant effect on the 


TABLE 5 
G*F Pole Placement and Robustness Design Pole Locations 





performance. The feedback matrix F is 
F = =0.27091 0.18869 0.0 0.0 
| 1.88298 -1.60017 0.0 0.00914 

The objective function is 0.0002 and the additive 
output singular value vs frequency is seen in Figure 6.21. 
This plot shows that in the robustness design the gain 
adjustment moves the minimum singular value from about 0.004 
rad, with very poor phase and gain margins, to a level of 
0.612 rad, close enough to the desired values of gain and 
phase. After obtaining the feedback yains the NPGS 'OPTSYS" 
program was used to obtain the necessary data, Bode plot, 
pole-zero map, and time response plot. The G*F robustness 
design time response plot is shown in Figure 6.22. The 
response of pitch rate is relatively good but sluggish 
response of pitch attitude to de step input. The performance 


of the pitch attitude was dejraded. It is possible that a 
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Figure 6.21 kobustness Design G*F Singular Value Plots. 


further change in the pole location couid achieve better 
performance and still meet the robustness criteria. The 
body axis forward velocity and vertical velocity perturba- 
tion responses are the Same aS in the pole placement only 
design(see Figure 6.4). Comparing the loop Bode plots of 
the pcle placement and robustness design transfer function, 
G*F, shows an increase in robustness. 

Considering a SISO Bode analysis, the channel input 1 
output 1 bode plot are shown in Figure 6.23 and 6.24. The 
bandwidth(BW) shows a small increase, 0.6 rad/sec, and the 
phase margin decrease of 26 degrees. The DC gain drops from 
21.4 db to 17.3 db, a4 db decrease and the slope change is 
Similar to the pole only design as seen in Figures 6.23 and 


6.24. This channel only shows a smooth curve. 
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Figure 6.22 G*F Robustness Design Time Response. 
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Figure 6.23 Pole Placement Only Design G*F 1:1. 
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Figure 6.24 Robustness Design  G*F 1:1. 
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Considering the channel input 1 output 2 with bode plots 
shown in Figures 6.25 and 6.26, a shift in the bandwidth(BW) 
from 2.4 rad/sec to 3.2 rad/sec is shown. The DC gaig 
increases from 12.7 db to 15.8 db at the minimum singular 
frequency of 1 rad/sec and an increase in the phase margin 
of 28 degrees. Figures 6.25 and 6.26 show the similar 
frequency response curve, Slope changes by zeros move the 
one pole position 0.8 rad/sec for pole-zero cancel and the 
other one Zero move to close the minimum singular value 
frequency 1 rad/sec and two pole location. The increase in 
the DC gain and phase margin yields an increased tolerance 
to perturbation. The channel input 2 output 1 bode plots 
are shown in Figures 6.27 and 6.28. The DC gain decreases 
16 db and the bandwidth shifts left 0.8 rad/sec, increasing 
the phase margin 151 degrees. The Zero shift has the effect 
oz smoothing the frequency response curve in the vicinity of 
the frequency of the minimum singular value. The channel 
input 2 output 2 bode plots are shown in Figures 6.29 and 
0, 308 These figures show a small shift left in the kand- 
width from 4.5 rad/sec to 2.8 rad/sec. The DC gain decreases 
12 db with an increased phase margin of 56 degrees. This 
results in a smooth curve, caused by the zero location shift 
in the same manner as previous channels(charaels 
1; 1,2; 1,2:2)- The G*F transfer functions Sole plot nut er- 
ical results are provided in Tatle 6. 

The noticeable change in the overall system, however, is 
in the transfer function input I Output 2. This is the 
change at the minimum singular vaiue position, increasing 
the gain and phase margin through a change in feedback gain. 
The optimizer routine brings the system gains to a more 
balanced condition and recovers a highly robust design. 

The gain changes associated with the robustness improve- 
ment cause the zeros of the various closed-loop pole-zero 


diagram of the closed-loop transfer matrix to move. A 
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Figure 6.26 Robustness Design  G*F 1:2. 
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Figure 6.28 Robustness Design  G*F 2:1. 
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Figure 6.30 Robustness Design  G*F 2:2. 
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TABLE 6 


G*F Transfer Functions Bode Plot Results 


channel Bandwidth eM 
E degree 
pole pole SV 





comparison of the eight  pole-zero diajrams is shown in 
Figures 6.31 to 6.34. 

The significant feature of these pole-zero diagrams is 
the shift of the zeros of the optimized design ina direc- 
tion that attempts to.equalize or balance the frequency 
response for frequencies in the vicinity of the minimun 
singular values. The pole-zero diagram of channel input 1 
output 2 will be discussed as an example of this effect. In 
figure 6.30, the pcle only design zeros are located about 
-0.31 and -0.76. When the pole placement and robustness 
routine completes the feedback gain modification these zeros 
EnXfr close to -1.7 and 0.28. The effect of these zeros 
shifting is to combine with the pole locations to equalize 
the frequency response and increase the DC gain 25 db at the 
same bandwidth and the minimum singular value frequency 1 
rad/sec, as depicted in Figures 6.35 and 6.36. Zero shifts 
for the remainder of the transfer functions provide similar 
results in the other channels. By moving toward the 
frequencies associated with the minimum singular values the 
zeros have balanced the overall frequency response of the 
system in each channel. While the channel gain modification 
is the primary mechanism for robustness recovery the zero 
shift associated with the feedback gain changes is directly 


related to the overall frequency response of the system. 
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4. 


Pole-Zero Map for 1:3 and 1 


Figure 6.32 
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2. 


Pole-Zero Map for 2:1 and 2 


Figure 6.33 





DESIGN 


ROBUSTNESS 


POLE ONLY DESIGN 


22 


4. 


Pole-Zero Map for 2:3 and 2 


Figure 6.34 
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Figure 6.35 Pole Placement Only Design G*F Closedloop 1:2. 
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TO summarize this analysis, a given performance level 
has been chosen in terms of pole locations. The level of 
robustness has been set for a desired gain and phase mar jin 
based on the universal gain and phase margin curve. The pole 
placement and robustness routine inproves the robustness 
level by changing the feedback gains that affect the channel 
input Í output 2 626055-03 ٣٦ This robustness recovery 
is affected by modification of the system feedback gains in 
such a manner that cross coupling gains are reduced so that 
small cross-coupling perturbations do not drive the systen 
iuto instability. The open-loop transfer function plots are 
used to indicate how this mechanism operates and have been 
Shown to be an alternative indicator of channels that may be 
affected by crossfeed perturbations. The pole-zero diajrans 
of the closed-loop transfer functions of the transfer matrix 
further indicate that zero movement is in a direction that 
equalizes the gain levei of the Frequency response curves in 
the vicinity of the lowest singular values providing a more 
balanced system response. ۶ ۰ a comparison of the 
input transfer function F*S output singular value and output 


transfer function G*F input singular value plots are shown 


in Figures 6.37 and 6.38. The F*G output ninimum singular 
value is too 10) 0.0076 rad). The G*F input minimum 
singular value is also its own low (0.1 rad). Therefore 


both transfer functions control only own function's robust- 
ness. Further analysis would be needed if robustness were 


required in both the input and cutput case simultaneously. 
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Figure 6.37 Robust Design F*G Output Singular Value Plots. 
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Figure 6.38 


VII. CONCIUSIONS 

The pole placement and robustness design routine coupled 
with the Automated Design Synthesis program provides the 
designer an excellent tool with which to attack the robust 
design problem. The pole placement and robustness design 
routine has demonstrated the capability of providing designs 
that solve the X22A V/STOL A/C longitudinal dynamic probiem 
caused by cross-coupling perturbations which reduce robust- 
ness in multivariable systems. This design improvenent is 
accomplished by modifying the system feedback gains in such 
a manner that the gain in channels that are affected by 
cross-coupling perturbations is eyualized with other systen 
gains to reduce this cross-coupling effect. The gain changes 
are accompanied by zero shifts which also influence the gain 
distribution and frequency respcnse of the system. 

The lack of robustness can be discovered in two ways. 
The first method is to plot the open-loop Bode plots of each 
element of the transfer matrix and look for extremely high 
gains or bandwidths relative to the other transfer func- 
tions. The second method examines the singular values of 
the return difference matrix for magnitude. The joint anal- 
ysis method tells where the robustness problem occurs in the 
design. Low singuiar values correspond to low robustness. 
The pole placement and robustness design routine can 
increase robustness by modifying feedback gains to reduce 
the effect of cross-coupling within the system. Observing 
the gain modification made by the pole placement and robust- 
ness routine the critical channel within the system that 
affects the robustness can be determined from the Bode 
plots. The pole placement and robustness routine feedback 


gain changes also cause zero shifts during the robustness 
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recovery. The gain on the open loop Bode plot for the 
affected cross-coupling channel is adjusted and the closes 
loop zeros, as seen on the pole-zero diayram, are Shitted. 
This zero shift is ina direction which will combine with 
System poles to smooth the frequency response liagram in the 
vicinity of the minimum singuiar value. The performance of 
the robustness design input transfer function F*y provides 
the best response, but the output transfer function  G*F 
indicates the performance of the pitch attitude is degraded. 
It is possible that a further change in the pole location 
would provide better performance and stili meet the robust- 


ness criteria. 
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